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Abstract 

We study the probabilistic evolution of a birth and death con- 
tinuous time measure-valued process with mutations and ecological 
interactions. The individuals are characterized by (phenotypic) traits 
that take values in a compact metric space. Each individual can die 
or generate a new individual. The birth and death rates may depend 
on the environment through the action of the whole population. The 
offspring can have the same trait or can mutate to a randomly dis- 
tributed trait. We assume that the population will be extinct almost 
surely. Our goal is the study, in this infinite dimensional framework, 
of quasi-stationary distributions when the process is conditioned on 
non-extinction. We firstly show in this general setting, the existence 
of quasi-stationary distributions. This result is based on an abstract 
theorem proving the existence of finite eigenmeasures for some pos- 
itive operators. We then consider a population with constant birth 
and death rates per individual and prove that there exists a unique 
quasi-stationary distribution with maximal exponential decay rate. 
The proof of uniqueness is based on an absolute continuity property 
with respect to a reference measure. 

Key words, quasi-stationary distribution, birth-death process, population 
dynamics, measured valued markov processes. 
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1 Introduction and main results 



1.1 Introduction 

We consider a general discrete model describing a structured population with 
a microscopic individual-based and stochastic point of view. The dynamics 
takes into account all reproduction and death events. Each individual is 
characterized by an heritable quantitative parameter, usually called trait, 
which can for example be the expression of its genotype or phcnotype. During 
the reproduction process, mutations of the trait can occur, implying some 
variability in the trait space. Moreover, the individuals can die. In the 
general model, the individual reproduction and death rates, as well as the 
mutation distribution, depend on the trait of the individual and on the whole 
population. In particular, cooperation or competition between individuals in 
this population are taken into account. 

In our model the set of traits T is a compact metric space with metric d. For 
convenience we assume diameter(T) = 1. Let B{T) be the class of Borel sets 
in T. The structured population is described by a finite point measure on T. 
Thus, the state space, denoted by A, is the set of all finite point measures 
which is contained in Ai{T), the set of positive measures on T. 

A configuration 77 G ^ is described by {riy : y G T) with rjy G Z_|_ = {0, 1, ...}, 
where only a finite subset of elements y G T satisfy rjy > 0. The finite set of 
present traits (i.e. traits of alive individuals) is denoted by 

{77} {y G T : ?7, > 0} 

and called the support of rj. For a function / defined on the trait space T, 
we will denote the integral of / with respect to 77 by 

ivJ) = f^y)^y 

Let I ■ I be the cardinal number of a set. We denote by ^7] = the number 
of active traits and by ||?7|| = ^ r]y the total number of individuals in rj. 

The void configuration is denoted by 77 = 0, so #0 = ||0|| = and we define 
:= ^\ {0} the set of nonempty configurations. 

The structured population dynamics is given by an individual-based model, 
taking into account each (clonal or mutation) birth and death events. 

The clonal birth rate, the mutation birth rate and the death rate of an individ- 
ual with trait y and a population 77 G ^, are denoted respectively by by{r)), 
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my{ri) and \{ti). The total reproduction rate for an individual with trait 
y G {q} is equal to hyirj) + myijj). We assume Aj^(O) = hy{^) = my{0) = 
for all y E T, which is natural for population dynamics. In what follows we 
assume that the functions 

Xyit]), by{r]), my(r]) : T x A^'^ R+ are continuous and strictly positive. (1) 

Let cr be a fixed non-atomic probability measure on {T,B{T)). The density 
location function of the mutations is : T x T :— IR+, {y, z) gy{z), where 
gy{-) is the probability density of the trait of the new mutated individual 
born from y. It satisfies 

J gy{z)da{z) = 1 for all ?/ G T . (2) 

T 

We assume that the function g.{-) is jointly continuous. To simplify notations 
we express the mutation part using location kernel G(?7, z) : B{T) IR+ 
given by 

Mr] e AMz e T , G{r],z) = ^ r]ymy{r])gy{z) = {ri,m.{v)9-{z)) . (3) 
Note that the ratio G(?7, z)da{z) / J G{ri, z)da{z) is the probability that, given 

T 

there is a mutation from t], the new trait is located at z. Hypothesis ([T]) and 
Lemma 11.41 stated below imply that the function G is continuous on ^ x T. 

We define a continuous time pure jump Markov process Y = (Yt) taking 
values on A. We denote by Q : A x B{A) M+, {r],B) Q{r],B), the 
kernel of measure jump rates given by 

Qiv,B)= Yl Yl J G{r],z)da{z). (4) 

The total mass Q{ri) of the kernel at rj is always finite and given by 

Q{'n)= Y{Q{'n.v+5y) + Q{'n.v-5y)) + jQ{v,v+Sz)daiz) . (5) 

y&{v} T\{,,} 
Observe that because of ([T]) 

VA; > 1 , Q+{k) = sup{g(r/) : \\ri\\ < k} < oo . (6) 
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The construction of a process Y with cadlag trajectories associated with the 
kernel Q, is the canonical one. Assume that the process starts from Yq = rj. 
Then, after an exponential time of parameter Q{r]), the process jumps to 
ri + 6y for y E {rj} with probability Q{rj,rj + 6y)/Q{ri), or to ri — 6y for y G {r]} 
with probability Q{ri, rj — Sy)/Q{ri), or to a point i] + 6z ior z E T \ {rj} with 
probability density Q{ri, t] + 5z)/Q{ri) with respect to a. The process restarts 
independently at the new configuration. 

The process Y can have explosions. To avoid this phenomenon and other 
reasons, throughout the paper we shall assume that 

B* = sup sup {by{r]) + my{r])) < oo. (7) 

This condition also guarantees the existence of the process {Yt : t > 0) as the 
unique solution of a stochastic differential equation driven by Poisson point 
measures. This is done in Section [2] following [H] , [5] . 

Since Q{0) = 0, the void configuration is an absorbing state for the process 
Y. We denote by 

To = mi{t > : = 0} 

the extinction time. In what follows, we will assume that the process a.s. 
extincts when starting from any initial configuration: 

Vr/G^: P^(To < cx)) = 1 . (8) 

So, in our setting we assume that competition between individuals, often due 
to the sharing of limited amount of resources, yields the discrete population 
to extinction with probability 1. Nevertheless, the extinction time Tq can 
be very large compared to the typical life time of individuals, and for some 
species one can observe fluctuations of the population size for large amounts 
of time before extinction ([TT])- To capture this phenomenon, we work with 
the notion of quasi-stationary measure, that is the class of probability mea- 
sures that are invariant under the conditioning to non-extinction. This no- 
tion has been extensively studied since the pioneering work of Yaglom for 
the branching process in [22| and the classification of killed processes intro- 
duced by Vere- Jones in . The description of quasi stationary distributions 
(q.s.d. for short) for finite state Markov chains was done in [U]. For countable 
Markov chains the infinitesimal description of q.s.d. on countable spaces was 
studied in [16] and [20] among others, and the more general existence result 
in the countable case was shown in [10] . For one-dimensional diffusions there 
is the pioneering work of Mandl [13] further developed in [3], [13], [IE] and 
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for bounded regions one can see [To] among others. For models of population 
dynamics and demography see [2], [12] and [3]. 

Let us recall the definition of a quasi-stationary distribution (q.s.d.). 

Definition 1.1. A probability measure v supported by the set of nonempty 
configurations is said to be a q.s.d. if 

V 5 G : F^iYt eB\To>t) = v{B) , (9) 

where B{A^^) is the class of Borel sets of A^'^ and where as usual we put 

When starting from a q.s.d. z/, the absorption at the state is exponentially 
distributed (for instance see ^0\)- Indeed, by the Markov property, the q.s.d. 
equality ¥^{Yt E dr],To > t) = iy{dr])F^{To > t) gives 

P,(To>t+s) = j ¥,{Ytedri,To>t+s) =¥,{To>t) j v{dr])¥^{TQ> s) 

_4-o _4-o 
= P,(To>t)P,(To>s). 

Hence there exists 0{v) > 0, the exponential decay rate (of absorption), such 
that 

Vt >0 : P^(To >t) = e-^('')*. (10) 

In nontrivial situations as ours, < Pv(7o > t) < 1 (for t > 0), then 
< eiu) < oo. 



1.2 The main results 

Let us introduce the global quantity 

A* = inf ini \y{ri). (11) 
»?e.4-o y&{v} 

Theorem 1.2. Under the assumption 

B* < K (12) 
there exists a q.s.d u, with exponential decay rate 

J Ml dnm 



This result is shown in Section |H It is based on an intermediate abstract the- 
orem proving the existence of finite eigenmeasures for some positive operators 
(Theorem 14.21) . 

In Section [5l we will introduce a natural a-finite measure fi and show that 
absolute continuity with respect to fi is preserved by the process. We study 
the Lebesgue decomposition of a q.s.d. with respect to fi. 

In Section [6] we will study the uniform case, which is given by 

Xy{r]) = A, by{r]) = 6(1 - p), my{r]) = bp, (13) 

where A, b and p are positive numbers with p < 1. The property f|T2|) reads 
X > b. In this case it can be shown that P = e~^^~''\ so Theorem 11.21 ensures 
the existence of a q.s.d. with exponential decay rate X — b. We will prove 
that this q.s.d. is the unique one with this decay rate, under the (recurrence) 
condition 

a0a{{y,z)eT^:gyiz) = O} = O, (14) 

and that given the weights of the configuration, the locations of the traits 
under this q.s.d. are absolutely continuous with respect to a. 

Theorem 1.3. In the uniform case assume that X > b and f|T4l) . Then 
there is a unique q.s.d. v on , associated with the exponential decay rate 
6 = X — b. Moreover v satisfies the absolutely continuous property, 

In this statement, ff denotes the ordered sequence of the elements of the 
support {rj}, (the compact metric space (T, d) being ordered in a measurable 
way, see Subsection 2.1), and 

f] = {r]y:y e {rj}) (15) 

is the associated sequence of strictly positive weights ordered accordingly. 

In all what follows, the set A will be endowed with the Prohorov metric 
which makes it a Polish space (complete separable metric space). This metric 
induces the weak convergence topology for which A is closed in the finite 
positive measure set. (See for example f7] Chapter 7 and Appendix). 

Let us give a general smoothness result which will be used several times later 
on. 
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Lemma 1.4. Let F : A x T —>■ M. be a continuous function on T. Then the 
function F defined on A x T by 

F{v, = j F{y,V,z)r]{dy) , 

T 

is continuous. 

Proof. Let ri,fj E A and z, z' G T. Thus 

\F{7],z)-Fifi,z')\ < {r],\F{.,r],z)-F{.,7],z')\) + {r],\F{.,7],z')-Fi.,fi,z')\) 

+ 1 {V -V,F{;V,z')) |. 

Since T is a compact set, it is immediate that the two first terms are small 
if z is close to z' and rj close to fj. If fj is in a small enough neighborhood of 
rj, these two atomic measures have the same weights, and the corresponding 
traits are close. In particular, fj belongs to a compact set, and the smallness 
of the last term follows by the equicontinuity of F on compact sets. □ 

2 Poisson construction, martingale and Feller 
properties 

Recall that ([H) and are assumed. We now give a pathwise construction of 
the process Y. As a preliminary result, we introduce an equivalent represen- 
tation of the finite point measures as a finite sequence of ordered elements. 

2.1 Representation of the finite point measures 

Since (T, d) is a compact metric space there exists a countable basis of open 
sets (Ui : i E N = {1, 2, ..}), that we fix once for all. The representation 

7^ : T ^ {0, 1}^ , z^ n{z) = (ci : i e N) with a = l{z e U^) 

is an injective measurable mapping, where the set {0, 1}^ is endowed with 
the product a— field. On {0, 1}^ we consider the lexicographical order 
which induces the following order on T: z < z' -vv- 'R.{z) <i TZ{z'). This 
order relation is measurable. 

The support {q} of a configuration can be ordered by ^ and represented by 
the tuple f] = {yi, ...,?/#^) and its discrete structure is = (ji{k) := rjy^ : k & 

k 

{1, #//}). Let us define So{r]) = and Sk{i]) = Yl Vyi for k G {1, #7]}. 

1=1 



7 



Remark that S^rjiv) = ll^ll- It is convenient to add an extra topologically 
isolated point d to T. Now we can introduce the functions if* : TU {d} 
by H^{ri) = for all rj & A and for i > 1 

H\7]) = ^ ^ (Sk-iiv), Skiri)] for A; < #r/ 

1 9 otherwise . 

The functions if* are measurable. We extend the functions b, X and m to d 
by putting bd{ri) = \d{ri) = mQirf) = for all r] E A. 



2.2 Pathwise Poisson construction 

Let {Q, JF, P) be a probability space in which there are defined two indepen- 
dent Poisson point measures: 

• (m) Ml {ds, di, dz, dO) is a Poisson point measure on [0, oo) x N x T x M+, 
with intensity measure ds {^^^^5k{di)) da{z)d6 (the birth Poisson 
measure) . 

• {i) M2{ds, di, dO) is a Poisson point measures on [0, oo) x N x M+, with 
the same intensity measure ds (^f^^^6k{di)) d9 (the death Poisson 
measure) . 

We denote {J-'t '■ t > 0) the canonical filtration generated by these processes. 

We define the process (Yj : t > 0) as a {J-'t : t > 0)-adapted stochastic process 
such that a.s. and for all t > 0, 



■t — -To 

[0,t]xNxTxn 



+ ^^l|f, a (z)<e<b - a (z)+m - a ,.-\\M^{ds,d%,dz,dQ) 

- / l{i<r.-ll}l{e<A^.(y^ )(y._)}^2(c?s,c/«,c/^). (16) 

[0,t]xNxR+ 

The existence of such process is proved in [TT], as well as its uniqueness in 
law. Its jump rates are those given by (jl]) so this process has the same law 
as the process introduced in Subsection II. 1[ In particular, the law of Y does 
not depend on the choice of the functions if* neither on the order defined in 
Subsection 12.11 
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Proposition 2.1. For any rj E A, any p > 1 and > 0, there exists two 
positive constants Cp and bp such that 

Er,{ sup \\Yt\\P) < Cp e^''*" < oo. (17) 

Proof. Let us introduce the following hitting times, 

V/ST e N : Tk = mf{t > : \\Yt\\ > K} . (18) 

From ffTUl) and neglecting the non-positive term, we easily obtain that for 
any K, 

WYtATj" < \\Yor+j{m-\\ + ^Y-\\ys-r)Mi{ds,dz,dz,de), 

Do 

where Dq is the subset of [0, t A Tk] x N x T x R.,. which satisfies i < \\Ys_\\ 
and 9 < bffi^Y,^)gH^{Y,^){z) + mHHYs^)9H-{Y,^){z). 

Then, by taking expectations, from ([71) and convexity inequality, we obtain 
E,(sup||yi^T,-r) < hr + 5>2^'-2 r(l + E,( sup \\Y^r))ds. 

t<to Jo u<sATk 

Standard arguments (Gronwall's lemma ) allow us to control the growth of 
the r.h.s. with respect to to- In particular, we have 

E^(sup \\Yt\\P) < {\\r]\\P + p2P'^B*)eP^'"'''*'\ 

With p = 1, (HI]) implies 

irP,(Tx <to) < hll e^**" , (19) 

and thus — *■ +oo a.s. as K —>■ +oo and the process is well defined on M+. 
Letting K go to infinity leads to the conclusion of the proof. □ 

Observe that the process is dominated everywhere by the integer- valued 
process Z solution of 

Zt=\\Yo\\+ J l{,<||z._||} 1|,<^. ^^^^^ ^^Mi(ds,ciz,c/2;,d^), (20) 

[0,t]xNxTxR+ 

which is a birth process with rate B*. This means that a.s. \\Yt\\ < Zt. We 
can establish the following stronger result. 
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Lemma 2.2. The process \\Y\\ is dominated by a birth and death process 
with birth rate B* and death rate X^. Then if we assume that \^ > B* , the 
process Y is absorbed exponentially fast. 

Proof. We introduce a coupling on the subset ^ oi Axfi defined by 

= {(r/,m) G ^ X N : \\r]\\ < m}. 

The coupled process is defined by its infinitesimal generator J, given by the 
rates 

J{ri,m;r] + 6y,m + l) = Vybyiv), V ^ {v} , 
J{r],m;r] + 6;,,m + l) = G{r],z), z ^ {r]] , 

J{r],m;r],m + 1) = mB* ~ J2 VyiKiv) + f^yiv)) , 

y^{v} 

J{r],m;r]-6y,m-l) = X^r]y , y e {i]} , 

J{r],m;r] - 6y,m) = Vy{>^yiv) V ^ {v} ^ 
J{ri,m;ri,m - 1) = K{m- Vy) ■ 

y<^{'n} 

It is immediate to check that the coordinates of this process have respectively 
the law of Y and the law of a birth and death process with birth rate B* and 
death rate A*. On the other hand when the coupled process starts from ^ 
it remains in ^ forever, so the domination follows. 

In [20] it is shown that the condition A* > B* implies that the birth and 
death chain is exponentially absorbed. The above domination implies that 
so does \\Y\\. 

□ 



It is useful to prove at this stage the following result on hitting times that 
only requires the property B* < oc. 

Lemma 2.3. For any t > and any rj G , there is a number c = 
ll^ll) ^ (0, 1) such that, 

Vir > : Fr,{TK<t) < c-^e-"^ . (21) 

Proof. The proof follows immediately from the domination of by the 
birth process Z introduced in ( 120|) . Indeed, assume Hloll < K and denote by 
T^j the smallest time such that Zt > M. Then T^ < T^ a.s.. Therefore, for 
any t > 0, and any t] G A"^ 

P,(Ti, <t) <P|Hi(r| <t) . 
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For a pure birth process (see for example [9]) we have 

%ii(r|<^)<PNi(^.>^) = EfrfA) • 

m=K ^ " ^ 

The result follows at once from this estimate. □ 
2.3 Martingale properties 

The process Y is Markovian and we describe its infinitesimal generator, in 
a weak form, using related martingales. The main hypotheses here are the 
boundedness of the total birth rate per individual (see ([7])) and the following 
bound for the death individual rate: there exist p > 1 and c > such that 

s\x^\y{ri) <c\\r]Y . (22) 

j/GT 

We define the weak generator of Y . Given / : ^ — M, a measurable and 
locally bounded function with /(O) = 0, we define Lf as 

Lf{r,) = Y,Vybyiv)UXv + Sy)-f{v)) (23) 

+ Vymyiv) j ifiv + Sz)-fiv))9y{z)da{z) 
ye{v} J 

+ E VyMv)ifiv-Sy)-fiv))- 

y&iv} 

Proposition 2.4. Let f : x A ^ be a measurable function such that 
for any p & A the marginal function f{*,p) is C^. We assume /(•,0) = 
and we take Yq = rj. 

(i) If f and dgf are bounded on [0,to] x A, for any to > 0, then 

J^l =: /(t, Y,) - /(O, r/) - [\dj{s, Y.,) + Lf{Y,))ds (24) 

Jo 

is a cddldg (JF^ : t > 0) -martingale, 
(a) Moreover, if there exists a finite p such that for any to > we have 

sup \f{t,v)\ + \dJ{t,7^)\<C{to){l+Mn, 
Q<t<to 

for some finite C{to), then is a martingale. 
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(Hi) If the functions f,dsf are assumed to be continuous, or more gener- 
ally locally bounded, then is a local martingale and for any T/v = 
inf{i: > : \\^t\ > ^} the process (-^/^/x^ '■ t > 0) is a martingale. 

Proof. Let us prove the first part of tlie Proposition. For alH > 

f{t,Yt)-f{o,r])- [ 9j(s,F,)ds = V(/(s,n_ + (n-y;_))-/(s,y;„)) 

"^0 s<t 

liolds almost surely. A simple computation shows that 

f{t,Yt)-f{0,v)- f dsf{s,Ys)ds = 
Jo 

[0,t]xNxTxR+ 

+ + l{e<^^^^^^^^iY,_)g^^^^^ ^)}^Mi{ds,di,dz,de) 

+ j{f{s, ^s-- 5H^{Ys.))-fis, Ys-))l{i<iiY,_ie<x^,^^^_^{Y,_)}M2ids, di, d9), 

[0,t]xNxK+ 

where both integrals belong to L^(P^). Compensating each Poisson measure, 
using Fubini's Theorem, and the fact that J gy{z)da{z) = 1, we obtain 

T 

f{t,Yt)-f{0,r])- [ {dsf{s,Ys) + Lf{Ys))ds 
Jo 

is a martingale. The rest of the Proposition is proved by localization argu- 
ments, justified by the result and the proof of Proposition 12.11 □ 

2.4 Feller property of the semi-group 

Let J\f = {Aft : t > 0) be the number of jumps for the process Y. We shall 
prove by induction the following result. 

Lemma 2.5. Assume that f : R_|_ xA^Wisa bounded continuous function. 
Then for all m > 

{t,ri)^E,if{t,Yt),N't = m) 
is a continuous function. 
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Proof. We notice that continuity and uniform continuity on every At, k > 1 
are equivalent because these sets are compact. Also we have that |/| is 
bounded on [0,to] x (Ufc=i -^fc) ^o,''^ and we denote by H/Hto,™ its 

supremum on this set. Denote hj i = \\ri\\ and n = m + i + 1. 

We first prove the continuity on time. For this purpose, we assume that 
< u < t < tQ where we assume that t, u are close and to is fixed. From 

/(t, Y,) = fit, F„)U=^„ + fit, Yt)l^,^j^^ 

we find (recall the notation (E])), 

\Kr,ifit,Yt),Aft = m) - Er,ifiu,Y^),Af^ = m)\ 

< sup \fit,^)-fiu,0\ + \\f\L,nQ+im + i)iit-u) + oit-u)). 

U\\<e+Tn 

Since the set {C|||C|| < i + is compact, it follows from the uniform 
continuity of / on compact sets that the first term on the r.h.s. is small if 
|t — m| is small. Hence the result follows. 

So in what follows we consider that t = u and we prove continuity on 77. We 
will do it by induction on m. In the case m = we have E^(/(t, Yt), Aft = 
0) = fit, r])e~^^'^^^ which is clearly continuous on i]. Now we prove the induc- 
tion step, so we assume that the statement holds for m and all continuous 
functions /. We have 



E^(/(t, Yt),Aft = mH) = / iA^ir], m, t-s}^A2iv, m, t-s}^A^ir^, m, t-s))e-'^^'^> ds. 



where 

Aiiri,m,t-s) = '^r]ybyir])Eri+Syifit- s,Yt^s), Aft~s = fn) 

yer] 

A2ir],m,t-s) = ^r]y\yir])En-5yifit-s,Yt^s),Aft~s = m} 
A^iri,m,t-s) = j En+5,ifit - s,Yt-s),M-s = 'fn)Gir],z)aidz). 

T 

Using Lemma 12. 5^ condition ([T]) and Lemma 11.41 it is immediate that the 
functions Ai, Ai and ^3 are continuous in (t, r/). We conclude by the Domi- 
nated Convergence Theorem since / is bounded. □ 
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Proposition 2.6. Let f : x A ^ M be a bounded continuous function. 
Then 

(t,r/)^E,(/(t,Fi)) 
is a continuous bounded function. 

Proof. Using Lemma (1) and the proof of Lemma r2.3[ we obtain tliat for 
eacli 7] E A, t > 0, there exists a = a{t, \\ri\\) > sucli tliat for any positive 
integer K, 

^v{Tk <t) = Pr,(Ar> M) <P^(Z > llr^ll) 

= P||^|| (Tf <t) < a-^e-"^ . (25) 

Assume tliat rj' is closed to ?7, and consider m, t close and smaller than to 
fixed. Then 

|E,(/(t,F,))-E,'(/(«,>^«))l <2||/||Pm(^,, >M+||77||) + 
f: |E,(/(t,ro, M = m) -E,,(/(M,y;), AC = m)\. 

m=0 

The result follows by taking a large AT, and by applying the bound ([25]) and 
Lemma I2.5[ □ 



3 Quasi-stationary distributions 
3.1 The process killed at 0. 

Let us recall that the state is absorbing for the population process Y . 
We have moreover assumed in ([H]) that the population goes almost surely to 
extinction, that is P(To < oo) = 1. This is in particular true if A* > B* . Our 
aim is the study of existence and possibly uniqueness of a q.s.d. z/, which is 
a probability measure on A~^ satisfying "^yiYt E B \ Tq > t) = vi^B). Let us 
now give some prehminary results for quasi-stationary distributions (q.s.d.). 

Since by condition ffT^ . the process Y is almost surely but not immedi- 
ately absorbed, and since starting from a q.s.d. the absorption time is 
exponentially distributed (see (fTOl) ). then its exponential decay rate satisfies 
< 9{-u) < oo. Since is absorbing it holds ¥^{Yt e B) = F^{Yt E B,To> t) 
for B G i3(^~°). So, the q.s.d. equation can be written as, 

Vfi G i3(^-°) , u{B) = e^(")*P,(rt G B) . (26) 
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From the above relations we deduce that for all 9 < d{v), Ej,(e^"^") < oo. 
So, for all 9 < d{v), u—a-.e. in rj it holds: E^(e^"^'') < oo. Then, a necessary 
condition for the existence of a q.s.d. is exponential absorption at 0, that is 

3t] G 36 > 0, E^(e^^°) < oo. (27) 

Let {Pt : t > 0) be the semigroup of the process before killing at 0, acting on 
the set Cb{A^^) of real continuous bounded functions defined on A~'^: 

yv e A-', V/ e CiA-') : (PtfM = ^r,ifiYt),To > t) . 

Let us observe that for any continuous and bounded function h : A ^ M. and 
for any rj G we have 

E^(/i(Ft)) = E^{h{Yt),To >t) + /i(0)P,(To < t) . (28) 

In particular, if h{0) = 0, we get E^{h{Yt)) = E^(/i(Ft),To > t). 

We denote by P/ the action of the semigroup on ^{A~^), defined for any 
positive measurable function / and any v G ^{A~^) by 

P!v{f) = v{PJ). 

From relation ( l26l) we get that a probability measure i/ is a q.s.d. if and only 
if there exists 6 > such that for all t > 

uiPJ) = e-%{f) , 

holds for all positive measurable function /, or equivalently for all / G 
CbiA''^). Then is a q.s.d. with exponential decay rate 9 if and only if 
it verifies 

Vt > : P/z/ = e-^V. (29) 
3.2 Some properties of q.s.d. 

Let us show that the existence of a q.s.d. will be proved if for a fixed strictly 
positive time, the eigenmeasure equation (1^ is satisfied. In what follows we 
denote by ^{A~^) the set of probability measures on A~^. 

Lemma 3.1. Let v G 3^{A~^) and /3 > such that P^D = (3i). Then (3 <l 
and there exists v a q.s.d. with exponential decay rate 9 := — log/3 > 0. 
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Proof. From /? = i>Pi(^-°) = P£,(To > 1) < 1, we get /3 < I, so 9 :-- 
— log/3 > 0. We must show that there exists u e ^{A~^) such that P^v - 
e"^* V for all i > 0. Consider, 



1 

e^'Plvds . 







For t e (0, 1) we have 



P}iy = fe'^Pl,uds= f \''Pl,uds+ f e'^P^^uds 
Jo Jo Jl~t 

/I rl+t 
e^(^-'^Pli)du + J e^^^'-'^Plvdu 

f e'^P^u du + e-^* /* e'^'e'PlPlv du = e-^V . 
Jt Jo 



For t > 1 we write t — n + r with < r < 1 and n e N. We have 



Note that 



□ 



0(iy) = lim - (1 - PJTo > t)) = lim - . (30) 

In the next result we give an explicit expression for the exponential decay 
rate associated to a q.s.d. Wc will use the identification between y & T and 
the singleton configuration that gives unit weight to the trait y. 

Lemma 3.2. If v & ^{A~^) is a q.s.d. then its exponential decay rate 9{i') 
satisfies 

0{u)^ J Q{r],0)u{dr]) = jQ{y,0)du{y) = J \y{y)du{y). (31) 

Proof. Since is absorbing we get that for all fixed G A^^ the absorption 
probability P,j(To < t) is increasing in time t. Let us denote 

a2{t) = sup{P^(ro < t) : \\r}\\ = 2} . 

Obviously we have 02(5) < a2{t) when < s <t. We claim that 

sup{P^(To < t) : ||r;|| > 2} < a2{t) (32) 
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Indeed let T2 = inf{t > : \\Yt\\ = 2}. Since the process will be a.s. extinct 
for all rj ^ A with \\ri\\ > 2, we have P^(T2 < 00) = 1. From the Markov 
property and the monotonicity in time of 02 (t) we get that for allr] ^ A with 

hll >2, 



?:||€ll=2-^° 

< «2(t) V / ¥r,if2 = ds,Yf^=0 = a2it). 



Now let us show that 02 (t) = o(t), that is lim a2(t)/t = 0. 

Let ?7 G ^2 be a fixed initial configuration, thus ||?7|| = 2. We denote by 
the subset of trajectories such that the function {\\Yt\\ : t < To) is decreasing, 
that is at all the jumps of the trajectory, an individual dies. Remark that, 
in the complement set of A^, either at the first or at the second jump 
of the trajectory, the number of individuals increases. Therefore, from fl32l) . 
the Markov property and the monotonicity in time of Pr,(To < t, Aj,), we get 
that 

sup{P^(To < t, A\) : llr^ll = 2} < sup{P,(To < t, A^) : Ht^H = 2} . 

Let us now denote by r the time of the first jump of the process F, and yi, 
7/2 are the locations of the points in rj (they can be equal). We have 

P^(To < t. Ay) < P^(To < t. Ay, K = yi) + P^(To < t. Ay, Y, = y^) , 

and 

¥ri{To < t, Aj, y^- = Ui) = P(e^ + By. < t, both events are deaths) = / fi{s)ds , 

Jo 

where and ey. are two independent random variables exponentially dis- 
tributed with parameters Qir]) and Q{yi) respectively. Moreover, condition- 
ally to the fact that the two jump events occur before time t, the probability 
to obtain two death events is x ~q^^- have for i = 1 (a similar 

computation holds for i = 2), 



/i(^)= / A,,(ry)e"«(^)%,(2/i)e-«(^^)(^-")rf«<Q(yi)(l-e-«(^)^). 
By using the bounds in ([6]) we find, 

sup{P,(To<t,Ai) : ||r/||=2}<Q+(l) /"(l-e-«+(2)^)cis = Q+(l) o(t). 

^0 
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So 02 (t) = o(t) holds and from flHUj) we obtain, 

9(1^) = \im J j Py(To < t)dv{y) . 

T 

Similar arguments as those just developed allow to get Py(To < t) = Xy{y){l — 
Q~Q{y)i^ + A; 02(^)5 where is a positive constant. Then the result follows. □ 



4 Proof of the existence of q.s.d. 

In this section we give a proof of Theorem II. 2[ The proof is based upon a 
more general result, Theorem 14.21 which shows that for a class of positive 
linear operators defined in some Banach spaces, whose elements are real 
functions with domain in a Polish space, there exist finite eigenmeasures. 
We show Theorem 11.21 in Subsection 14.21 For this purpose we construct the 
appropriate Banach spaces and the operator, in order that the eigenmeasure 
given by Theorem 14. 21 is a q.s.d. of the original problem. 

4.1 An abstract result 

In this paragraph, {X, d) is a Polish metric space. We will denote by Cfe(A') 
the set of bounded continuous functions on X. This set becomes Banach 
space when equipped with the supremum norm. 

Let 5* be a bounded positive hnear operator on Cb{X). We will also make 
the following hypothesis. 

Hypothesis J^: There exists a continuous function on {X, d) such that 
M[ ^2 > 1. 

Jif2 For any u > 0, the set v'^"'^([0, u]) is compact. 

It follows from that if {X, d) is not compact, there is a sequence {xj : j G 
M) in X such that limj^oo 9^2 (a^j) = 00. 

Before stating the main result of this section we state and prove a lemma 
which will be useful later on. 

Lemma 4.1. Let v be a continuous nonnegative linear form on Cb{X). As- 
sume there is a positive number K such that for any function ip G Cb{X) 
satisfying < ip < (p2, we have 

v{ip) < K . 
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Then there exists a positive measure v on X such that for any function f e 

v{f) = j fdv. 

Proof. Let Cq{X) be the set of continuous functions vanishing at infinity. Let 

w be a real continuous non- increasing nonnegative function on M+. Assume 
that w = 1 on the interval [0, 1] and w{2) = (hence zu = on [2, oo)). 
For any integer m, let Vm be the continuous positive linear form defined on 
Co{X) by 

Vmif) = v{'nj{(p2/m) f) . 

This linear form has support in the set (P2^{[0, 2m]) in the sense that it van- 
ishes on those functions which vanish on this set. Note also that <^^^([0, 2m]) 
is compact by hypothesis J^2- Therefore it can be identified with a nonneg- 
ative measure on X, namely for any / e Ch{X) we have 

M/) = 

We now prove that this sequence of measures is tight. Let u > and define 
the set 

Again, by hypothesis J^, for any u > this is a compact set. We now 
observe that < 1 — zu{2(f2/u) . Therefore, 

^m{K^) < Vm{^ - W{2(p2/u)) = Vm{l " Zu{2(p2/u)) 

= v{w{if2/m) (l - 07(2 v?2/m))) , 

We now use the fact that the function ip2w{i~p2/m) (l — w{2ip2/u)) is in 
C}){X) and satisfies 

11 

-w{tp2lm) (l - w{2 ^p2/u)) < w{Lp2/m) (l - w{2 tp2/u))tp2 < ^2 

to obtain from the hypothesis of the lemma that 

2 2K 

v{w{Lp2/m) (1 - cc(2(^2/m))) < -v{uj{Lp2/m) (l - 'uj{2lp2/u))lp2) < . 

In other words, for any u > we have for any integer m 
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The sequence of measures Um is therefore tight, and we denote by u an 
accumulation point which is a nonnegative measure on X. We now prove 
that for any / G Cb{X) we have z/(/) = v{f). For this purpose, we write 

v{f) = v{zu{ip2/m) /) + v{{l - w{ip2/m)) f) . 

We now use the inequahty 

> (1 - ■C0{ip2/m))ip2 > "^(1 - , 

to conclude using the hypothesis of the lemma (since (1 — 'Dj{ip2/m))ip2 G 
Cb{X)) that 

\v{{l-w{^2/m)) f)\ < v{{l-w{^2lm)) |/|) < ||/|| v{l-w{^2/m)) < ^ . 
In other words, we have for any / G Cb{X) 

\vU)-^mU)\<- ■ 

' ' m 

From the tightness bound, we have for any / G Ch{X) 

lim z/„(/) = z/(/) , 

see for example [T], and therefore v{f ) = v{f) which completes the proof of 
the lemma. □ 

We now state the general result. 

Theorem 4.2. Assume hypotheses M'\ and M2, and assume also that there 
exist three constants Ci > 7 > and D > such that 

3(1) > ci 

and for any ip G Cb(A') with < ip ^ f2 

Sip < 7V92 + D . 

Then there is a probability measure u on X such that u o S = I3p, with 
13 = z/(5(l)) > 0. 

Proof. In the dual space Cb{X)*, we define for any real K > the convex 
set J(fK given by 



jeK= {ve Cb{X)* : v>0, v{l) = 1, sup v{tp) < K 
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Note that by Lemma 14.1^ the elements of are positive measures. 
We observe that for any K large enough the set is non empty. It suffices 
to consider a Dirac measure 6x on a point x ^ X and to take K > y:>2{x). 
Since for any K > 0, is an intersection of weak* closed subsets, it is 
closed in the weak* topology. 

We now introduce the non-linear operator T having domain and defined 
by 

V o S 

Note that since 5(1) > ci > 0, we have f (5(1)) > civ{l) and this operator 
T is well defined on .^k- We have obviously T(t>)(l) = 1. We now prove 
that T maps J^k into itself. Let ip G Cb{X) with < < ip2- Since 

Sip < 7V?2 + D 

and obviously 

0<Sip<-f 



7 

we get 

0<SiP<^{^2Ai\\SiP\\h))+D. 

Therefore since the function ip' = A {\\Sip\\/'y) satisfies ip' G Cb{X) and 
<■?/''< 9?2- We conclude that for v G 



Cl 

From the bound v{ip') < K we get 

T{v) {^) < ^ < ±K + — < K 

Cl Cl Cl 

if i^' > D/(ci — 7). Therefore, for any K large enough, the set is non 
empty and mapped into itself by T. 

It is easy to show that T is continuous on in the weak* topology. This 
follows at once from the continuity of the operator S. We can now apply 
Tychonov's fixed point theorem (see |I9| or |8j) to deduce that T has a fixed 
point. This implies that there is a point u G J(^k such that v o S = v{S{l))v. 
This concludes the proof of the Theorem. □ 
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4.2 Construction of the function Lp2 and the proof of 
Theorem 

We assume that the hypotheses of Theorem 11.21 hold. In our apphcation we 
have a semi-group Pt acting on Cb(^~°). We will use Lemma [HTT] to construct 
a q.s.d. This lemma is proved using Theorem 14.21 applied to 5* = Pi: 

Sfiv) = Pifiv) = , To > 1) , G 

Here the Polish metric space {X, d) of the previous paragraph will be {A~^ , dp), 
and so the function (f2 will have domain in the set of nonempty configura- 
tions. 

We recall the elementary formula valid for any continuous and bounded func- 
tion / on ^ and any t >0 

E,(/(Fi)) = E,(/(Fi) , To > t) + /(0)P,(ro < t) . 

We start with the following bounds. 

Lemma 4.3. Let Ai = sup sup Xy{ri) < oo, then 

T): \\r)\\=l y&l 

-Ai < LI < . 

and for all t > 0, 

g-Ait < p^i < 1 . 

Proof. The proof follows at once from Lemma 12.41 and a computation of 
Ll^-o (see formula (1^ ). □ 

Lemma 4.4. Consider for any a > the function ^2{''l) = e'^"''" l_4-o (r^) . 
Then 

L^2ir]) < {B* {e'^ - I) + K (e^ - l)) U\ v.^(r/). 

Proof. We compute Lip2{ri) using (1231) . For t] G we have 
Lifiiv) = E r/,(6,(r/)+m,(r7))(e»-l)e»W 

+ E riyXyiv) {e-^ - 1) e^M _ A,(r7)l||,||=i (33) 

y^{v} 

< B* \\r]\\ yD^{r]) (e" - 1) + A, \\r]\\ (/?^(r/) (e"*^ - 1) . 

□ 

To define the function ip2, we will need the two following results. 
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Lemma 4.5. The differential equation 

^ = A,(l-e-'^)+S*(l-e'^) (34) 

has two fixed points a = and a = log(A*/-B*). The trajectory of any 
initial condition ao G (0, log(A*/i?*)) is increasing in time and converges to 
log(A,/5*). 

Proof. Left to the reader. □ 

Lemma 4.6. Assume ffT^ . Let a{t) he the solution of with initial 
condition qq G (0, log(A^,/i?*)). Then 



sup E^(e-^**e"WII^*ll,To > t) < e 
teM+ 

Proof. We introduce the function 

/(t,r^) = e-^**e"Wll^lll^-o(r^), 

and for any integer N we denote by f^ the function 

r(t,r/) = /(t,r/)l|,,||<^. 

Note that f^{t.,ri) is continuous with compact support {rj : \\ri\\ < N}. 
Using Proposition 12.41 {iii) we get 

ftATM 



f%t,Y,^T^) = f {Q,Y,) + / 

^0 



where is a martingale. Then we obtain 

/ rthTM 

Observe that if > M and s < Tm we have = /(s.^^)- Let 

tend to infinity to get 

/ rt/\TM \ 
E,(/(t,WM)) =/(0,lo)+E, ^ {dsf{s,Y,) + Lf{s,Y,))ds \ . 

Using Lemma [4.41 we have 
dJ{s,Ys) + Lf{s,Y,) = 

e-^- ((A. (l-e-"(^))+i?* (l-e'^(^))) r.||-A.) vf\Ys) + Liff\Y,) < 0. 
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Therefore 

</(0,>o) ■ 

Letting M tend to infinity and by using the Monotone Convergence Theorem 
we obtain, 

</(0,Fo) . 

The result follows from the definition of / □ 
We take as function ip2 the function 

a(l) 

for a solution of ( l34l) with initial condition ao G (0, log(A*/i?*)). The operator 
S is given by S* = Pi, and hence is positive and maps continuously C;,(^^°) 
into itself. 

We must now show that S = Pi, and satisfy the hypothesis of Theorem 

Mi 

Lemma 4.7. 

(i) The hypotheses Jifi and M2 are satisfied. 

(a) S{1) > ci > 0, with ci = e"^i. 

(Hi) For any 7 > 0, there is a constant D = -D(7) > such that for any 
^ G Cb{X) with < ^ < 

Sil) < + D . 

Proof. The hypotheses and J^2 are easy to check using the Feller property 
of Pi (see Proposition 12.61) . 

{ii) follows at once from Lemma [4.31 We now prove {Hi). 
Let ip G Cb{X) with < < ip2. We have from Lemma [4.61 

PiiP{7^) = E,(^(ri) , To > 1) < E,(^2(yi),To > 1) < e"* e"""^" . 

Since a(l) > ao by Lemma [4.51 for any 7 > there is an integer such 
that for any m > we have 

Therefore, for any rj we have 

In other words, we have proved {Hi) with the constant D = e^* e°'°"^^. □ 
Theorem 11.21 follows immediately from the previous Lemma and Theorem 
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5 The process and absolute continuity 



In this section we introduce a natural a— finite measure We will show that 
the process Y preserves the absolutely continuity with respect to /i and that 
when the process starts from any point measure after any positive time the 
absolutely continuous part of the marginal distribution does not vanish. 



5.1 The measures 

We will denote by T'^ the set of all A;— tuples in T ordered by ^ defined in 
Subsection 12.11 So, for r] E A, its ordered support ff = ?/#,,) belongs to 

f #'?. The discrete structure f] is an element in N*'' and the set of all discrete 
structures is denoted by 

Here N° contains a unique element denoted by and it is the discrete struc- 
ture of the void configuration rj = 0. A generic element of E(N) will be 
denoted by q. Moreover for each q e we put #g = k if qEN^. We put 



A^= {r] e A : 1] = q} for g G S(N) 



and for B d A, 



= {t] e B : f] = q} for g G S(N) . 

In the sequel for g G and C C T'^ we denote 

{q} X C := {i] e A:r] = q,f] e C} . (35) 

We denote by ^f{A) the set of measures on {A, B{A)) that give finite weight 
to all sets Ak- By ^/(S(N)) we mean the set of measures on S(N) giving 
finite weight to all the subsets N'^, and ^f(N)) denotes the measures on N 
giving finite weight to all its points. Every measure v G ^f{A) defines a 
measure v G by 

v{q) = v{A^) . (36) 

Also V defines a set of conditional measures Vg- G ^ (¥#"?) by 

ifv{q)=0, 
v{r] E A:r] = q, f] G •)/v{q) otherwise . 
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Then e ^(T#9) is a probability measure if v{q) > 0. In the case that 
V e ^{A) we have 

M*) ^v(f]e»\f]^q) . (37) 

Conversely a probability measure G ^(^) is given by a probability measure 
e ^(E(N)) and the family of conditional measures {vjj e ^(^7^)) so that 

v{B)= mv.iB^), BeB{A-'). 

geS(N) 

In this sense 

dv{r]) =v{rf)dvrj{ri). (38) 

Let (/9 : ^ — > M be a function. Observe that its restriction to Afj can be 
identified with a function 9?! with domain in by the formula . {ff) — 
ip{r]). Let : ^ ^ R be a integrable function, we have 




Now we define the measure jji by, 

//(N^^ X TO) = /.({O}) = 1 and = x ^ for > 1 , (39) 

where is the point measure on N'^ that gives a unit mass to every point, 
and is the restriction to T*^ of product measure cr®^. Note that v e ^f{A) 
satisfies 

V « n ^ (yq^ E(N) : « a*?^ . 
Hence, if v e J^{A) is such that v « /i, then v is of the form 

v{r] e A -.f] = q,f] e dy) = v{q) (p^{y) da*'i{y) . (40) 

where v G ^(E(N)) and for each fixed q G E(N), ^q{») is a density function 
in T#9 with respect to (7#9. 

5.2 Absolutely continuity is preserved 

Proposition 5.1. The process Y preserves the absolutely continuity with 
respect to fi, that is 

yve ^(A), v«ii ^ > ¥^{Yt e •) « //(•) . (41) 
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Proof. Let us define the jump time sequence, 

To = and r„ = inf{t > : Yt ^ Yr,,-^ for n>\. 

In particular r = ti is tlie time of tlie first jump. Remark tliat tlie sequence 
r„ tends a.s. to infinity, as it can be deduced from ffTUl) . Wlien we need to 
empliasize tlie dependence on the initial condition Iq = ^ we will denote 
and instead of r„ and r, respectively. We have 

Since the sum of absolutely continuous measures is also absolutely continuous 
it suffices to prove that 

P^,(yi G •, r„, < t < r„,+i) « II for all n > . (42) 

First, let us show the case n = 0. For B G we have that the expression 

P,(rtG5,t<r)=P,(roGS,t<r)= / i;(dr/)P^(t < r) (43) 

vanishes if viB) = 0, so also when fi{B) = 0. 

Before considering the case n > 1 in fH21) let us prove the relation 

V « fi ^ P„(y, G •) << /i. (44) 

It suffices to fix g* G S(N) and to show that the measure P„(Fr £ •) restricted 
to the class of sets ({g} x C : C E i3(T#'?)) is absolutely continuous with 

respect to cr#'? (see [351). Let = For A; = the claim holds because 
/x({0}) = 1. Let k > 1. Recall that the probability measure v has the form 
stated in fHOl) . so iprj = dvjj /da'f^'^ is the density function on the space T#^. 
Below, for ?7 G ^ we denote 

r^~y = rj - 5y, y e {r]}; ?7+^ = r/ + 5^ , V2; G T ; 

r^-y rj+^:=^; fj-v ■= r]~y ; := 77+^ . (45) 

We denote v^ = cr^(T^) (for A; = we set vq = 1). Let A; G N, g G and 
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C e B(T^) be fixed. We have 

P.(n e {q} ^C) = J dvii) J -^^Q{r^\dv) 



\yG{'n}-r, 

+ /i(^=^)( E 

c \y^{v] 



{q}xC 

Qiri-y) 

(r/,-l)A,(r^+^) 



+ 



/ 

l{r] = q) 



c 



\zeT\{v} 



Qiv 



da^ij]) 



J 



\da''{f]) 



where in the last two terms we have used the following relations 

Vfc 



da{z)da''{f]) 



and 



da^{j]) 



Vfc-l 



Hence, the relation (l44l) is proved. An inductive argument gives 
P^(y^^^ G •) << for every n > 1 . 



(46) 



Now, let us show that (l42l) holds for n > 1. Denote by F„ the distribution 
of r„. By the strong Markov property and Fubini theorem we get 

^.{Yt e;Tr,<t< r„+i) = [ (Py^^ {Yt_, g •, r > t - s)) . 



which, by using relations ( H3l) and ( H6l) . is absolutely continuous with respect 
to fi. □ 



5.3 Evolution after the first mutation 

We want to study the absolute continuity with respect to /i of the law of 
Yt initially distributed according to a general measure v. To this aim we 
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will introduce the first mutation time. Note that a mutant individual has a 
different trait from those of its parent, so the time of first mutation is 

X = inf {t > : {Yt} ^ {Fq} } • (47) 

When X is finite we have {Y^} ^ 0, so (x < oo) =^ (x < Tq). 

Now, let us consider the first time where the traits of the initial configuration 
disappear, 

K = mi{t > : {Yt} n {Fo} = 0} , 

and for a fixed r], the first time where the traits of rj disappear, k'^ — inf{i > 
: {Yt} n {rj} = 0}. When {rj} n {Yq} = 0, then k1=0. 

We have k x except when k — x — oo. Obviously k <To. Moreover 

(x>k) <^ (oo = x>«) <^ (x>« = 7o) and {kkTq) ^ {x<k<To) . (48) 
Also note that (k < x) H (k < t) C = Tq < t). Since P^(Yf G •, Tq < t) = 
So{») is concentrated at = 0, then 

p^(rte.,«<x,«<i) = 5o(«). 

The unique nontrivial cases are the following two ones. 
Proposition 5.2. Let r] e and t >0, we have: 

(i) V'n{Yt E X < 1^ ^ t < Tq) is absolutely continuous with respect to fi and 
it is concentrated in A^^ ; 

{ii) P^(yt e < k) is singular with respect to /i. 

Proof. Let us show (i). Prom the Markov property we have, 

¥,{Yte;x<^i<t) = Yl ^vix<f^<t,Y^-=^,Yte*) (49) 

J gy{z)F^+.{Yt_, e •, K^<t-s)da{z). 

Hence, it is sufficient to show that for every u > 0, rj & A~'^ and y e {77}, it 
holds 

j P^+.(y; u)gy{z)da{z) « 
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By using f P^+z {{^u} H {rj} 7^ 0, k'' < u) gy{z)da{z) = 0, and since the 

measure a is non-atomic, a similar proof to the one showing Proposition 
15.11 works and proves the result. Indeed, for each t > 0, the singular part 
with respect to /i of fr^{Yt G ■) is a measure on the set of atomic measures 
with support contained in {rj} (corresponding to death or clonal events from 
individuals initially alive). 

Let us show (ii). Let {rj} C T be the finite set of initial traits and put 
k = #77. Consider the Borel set B = e : {^} n {r]} 7^ 0} and define 
Bi^n = e ■■M = n, |{0 n {v}\ = /} for n e N, Z = l,...,n A fc. We 
have B = |J -B/,™- Since a is non-atomic we have fi{Bi^n) = for 

neN, /e{l,...,nAfc} 

all / G {1, ...n A k}. On the other hand, from the definition of k we have 
Fr,{Yt e B,t < k) = 1, and the result follows. □ 

Let V G ^{A). We denote by f* the distribution of Yt when the distribution 
of Yq is V, that is 

v\B) = F^{Yt G 5) , Be B{A) , t > . (50) 

We denote by f = v^'^-i-v^^ the Lebesgue decomposition of v into its absolutely 
continuous part v^^ << fi and its singular part v^^ with respect to /i. For f* 
this decomposition is written as f * = f ^'^'^ + v^'^\ As usual, 5^ is the Dirac 
measure at 77 G ^, so 5* denotes the measure ^^(•) = ^riiXt G •). We will 
denote by supp(f ) the closed support of a measure v. 

Proposition 5.3. The process Y verifies: 

(i) For allt>0 and all rj G we have (5*'''''(^"°) > 0; 

(ii) For allt>0 and all v G ^(^-°) it holds v^^"" > J^.,, S^^^^vidr]) > 0; 

(in) Assume condition [14\ )- Then for all rj G A"^ with {77} C Supp{a) and 
for all e > 0, the following relation holds, 

Vt>0, 5;'"^(i?(r7,e))>0, 

where B{rj, e) = {77' G ^ : \\rj' — rj\\ < e}. 

Proof. It suffices to show {i) and {Hi). Let us show the first part. Fix t > 
and rj G A~^. We claim that P^(x < < ^ < To) > 0. In fact, it suffices 
to consider the event where a mutation occurs at the first jump and after it 
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all the initial traits disappear before t and these changes are the unique ones 
before t. This event has strictly positive probability, so the claim is proved. 
Proposition 15.21 (i) gives fniXt ^•,X<^^^< ^o) << and we deduce, 

(5;'^^(^-°)>P^(x</t<t<To)>0. 

Then (i) holds. 

The proof of {in) is entirely similar to the proof of {{) but we need some 
previous remarks. For every ?/ G T we have J gy{z)da{z) = 1, and so, 

T 

(j{{z G Supp(cr) : gy{z) > 0}) > 0. On the other hand from condition ( 1141) 
the set 

D = {ye Supp(a) : a{{z G Supp(a) : gy{z) > 0}) = 1} 

verifies <j{D) = 1. In particular cr(z G -D : > 0) > is satisfied for all 

?/ G T. Now, let {t]} = {ui : i = l,...,k} and consider the following event: 
a mutation occurs at the first jump to a trait y' G D, afterwords successive 
mutations to the traits in B{yi, e)nD take place, then for each trait yi there 
are Qi — l clonal births, and finally all the initial traits and y' disappear. This 
history occurs before t and assume that these changes are the unique ones 
that happen before t. 

From condition fll4p and the definition of D this event has strictly positive 
probability and the claim is proved. □ 

5.4 Decomposition of q.s.d. 

Let us study the Lebesgue decomposition of a q.s.d. with respect to /i. 

Proposition 5.4. Let u be a q.s.d. on A^^. Then, 
(i) u'^^ ^ 0; 

(m) Assume Condition [T^ . Then {rj G A^^ : {i]} C Supp{a)} C Supp{v°''^); 
{Hi) If v^^ 7^ 0, the probability measure v*'^^ := z/*Yz/'*'(^^*') satisfies 

F^*4Yt e B) = e-^'^''^'u*''{B) ^ B e B{B'') ,t > , (51) 
where 5*' G B{A^^) is a measurable set such that ii{B'^^) = and z/**(i?**) 
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Proof. We first note that the existence of the set B^^ is ensured by the Radon- 
Nikodym decomposition theorem. Set H := A^^ \ Let us show that 

yt>0,\/7] e A-^ : 6l{H)>0. (52) 

Since fi{B'') = and 5*'^^ « fx we have 5*-^ (5^^) = 0. Then 5l;^%H) = 
§tAc(^j^~oy gy Proposition Ea > for alH > and all t] e A'^. 

So 

6l,{H) > 5f\H) = > 0, 

and the assertion fl3^ holds. 

Now we prove part (z). We can assume ^ 0, if not the result is trivial. 
From fl^ we get, 

u\H) = [ 6l{HHdv) > 0. 

On the other hand, from relation (1261) we obtain z/(if) = e^^^^^u^^H) > 0. 
Since v^^{H) = 0, we necessarily have v^{H) = z/(if) > 0, so (i) holds. Now, 
from Proposition 15.31 [in) a similar proof as above shows {ii). 

Let us show (iii). Let u*'^ := u"^ {A'^) . For every B C B G B{A-°), 
we have 

z/(5) = e^^'^)* (z/'^=(^-°)P,.ac(yi e 5) + z/^X'4"°)P.*4^t e B)) . (53) 

By Proposition 15. Y preserves /i, so P,^*ac(Y'i g •) << /i. Since fi{B^^) = 
we get ¥^*.c{Yt G 5"^) = 0. By evaluating (p]) at t = and since B G i3(5''') 
we find v*^^{B) = v{B)/v{B^^). By putting all these elements together we 
obtain relation ( l5Ti) . □ 

6 The uniform case 
6.1 The model 

In this section, we assume that the individual jump rates satisfy, 

\y{ri) = A , by{ri) = 6(1 - p) , my{r]) = bp, \fy e {r]} , 

A, b and p are positive numbers with p < 1. Recall that g : T x T 
a jointly continuous nonnegative function satisfying J gy{c)da{c) = 

T 

y eT and the condition ( fT4|) . 
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We observe that in this case the process of the total number of individuals 
||y|| = {\\Yt\\ : t > 0) is a Markov process and that Yt = \\Yt\\ = 0, which 
means that the time of absorption at of the processes Y and is the 
same (note that even if the O's have a different meaning, they are identified). 

Now, in [20] it is shown that there exists a q.s.d. for the process killed 
at if and only if A > 6. In addition, the extremal exponential decay rate 
of defined by sup{^(i/); u q.s.d.}, is equal to A — 6 and there exists a 
unique (extremal) q.s.d. C for with this exponential decay rate X — b, 
given by 

C{k) =(J^ (^I-jY k> 1. (55) 

When z/ is a q.s.d. for Y with exponential decay rate ^(z/) then the probability 
vector ( = {({k) : k ^ N) given by 

Cik) = iy{Ak), keN, (56) 

is a q.s.d. with exponential decay rate 6 = Oiy), associated with the linear 
birth and death process ||V||. Hence a necessary condition for the existence 
of q.s.d. for the process F is A > 6. We also deduce that all quasi-stationary 
probability measures i) oiY with exponential decay rate A — 6 are such that 
v{Ak) = so by ([55]) we get 

z>(9?i) < oo, where ipi{ri) = \\ri\\. 

Now, we know from Theorem 1 1 . 21 that there exists a q.s.d. u with exponential 
decay rate ^(z/) = z/(Pi(l)). Moreover, it is immediate to show that (fi 
satisfies Lipi = —(A — b)ipi. Then, from Proposition 12.41 we get PiV^i = 
g-(A-b)^^_ Hence, if z^ is a q.s.d. provided by Theorem 11.21 its exponential 
decay rate should be 

e = x-b. 

Let us now consider the semi-group Rt given by, 

RtivM = e^'PtivM = e'%,{^{Yt)lT,>t) , t>0. 

The function (pi satisfies Rt<Pi = (pi u—a.e. for all q.s.d. z/ with the expo- 
nential decay rate 6. 

Proposition 6.1. Every q.s.d. v with exponential decay rate 6 = X — b is 
absolutely continuous with respect to fi. 
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Proof. Let z/ be a q.s.d. which is not absolutely continuous. Then we can 
write the Lebesgue decomposition 

j, = f^^^ thatisz/(5)= I fdfi + ^{B), B eB{E), 

Jb 

where / is a nonnegative /i— integrable function and ^ is a singular measure 
with respect to fi. 

From now on we denote by Rj the dual action of Rt on the set of measures 
defined by (i?Jw)((y9) = v{Rt^p) for every measure v G ^f{A) and any positive 
measurable function Lp. Since v is a. q.s.d. it is invariant by the adjoint semi- 
group i?J, that is = z/, then 

On the other hand, it follows from Proposition 15. II that R\{f n) <^ /x. There- 
fore 

Since (pi is u integrable it must also be fdfi integrable. From the relation 

Rti^pi) = fi we get, 

j y^i fdfi = J <^idRl{ffi), 

and since ipi is strictly positive, we conclude R\{f ^) = ffi. This implies 
RtiO = C However by Proposition 15.31 (ii) (with v = C,), -RJ(0 cannot be 
completely singular with respect to /i unless ^ vanishes. This concludes the 
proof of the proposition. □ 

Let us now turn to the study of uniqueness. 

Lemma 6.2. Assume condition [T^: cr (g) (T{{g = 0}) = 0. Then the Borel 
set ^(1,1) = {rj E A : q = {1, 1)} satisfies fi{A(i^i)) > 0. Moreover, for any 
q.s.d. V with exponential decay rate 6 = X — b and for any Borel set B with 
nB) > 0, we have z^(^{i,i) CiB) > 0. 

Proof. Let us consider a q.s.d. u with exponential decay rate 9 = \ — b. 
Lemma [32] implies that the restriction of to .Ai = {?7 G ^ : \\ri\\ = 1}, 
does not vanish. On the other hand by the previous result, it is absolutely 
continuous with respect to a. Then, 

di^(i) = fida 



34 



for some nonnegative function /i that does not vanish on a set of a positive 
measure. 

For any function / in Ci,{A) such that /(O) = and such that is supported 
in a compact set, that is /(?/) = for all ||?7|| large enough, it follows from 
v{Ptf) = exp(-^t) v{f) that 

u{LfU-o) = -e u{f) . 

Since this is true for any such function, we get (with rj = {rjy : y E {q}) as 
defined in ffT^ . and notation flHUl) and 



-edu^{f^) = h{l-p) {Vy~mV-')diyjj-y{v) 

y:r)y=l y' e{r)}\{y} 

-(A + 6)l J^r^J mdMv). (57) 
\y&{v} ) 

It follows from equation ( 1571) applied to the measure v and solving for that 
for some constant C > we have 

dy{i,i){yi,y2) > <^ (/i (1/1)^2/1(1/2) + /i (1/2) ^2/2(^/1)) da{yi)da{y2) . 

Using this lower bound in the equation for z/(i), we get for some constant 
C >0 



fM>C' j h{u)gMda{u). 
Therefore for some constant C" > we have the estimate 
dT^(i,i){yi,y2) > C" daijji) da{y2) x 

9yAy2) I fi{u)gu{yi) da{u) + gy^{yi) l fi{u)gu{y2)da{u] 



Let 5 be a Borel set with ^{A{i^i) fl i?) > 0. By the identification between 

^(1^1) and (1, 1) X it can be assumed that B C T^. We get from Fubini's 
Theorem 

1^(^(1,1) nfi) > C" j lB{yi,y2)gyAy2) fi{u) gu{yi)da{u) da{yi)da{y2) . 
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Therefore, if ^{A^i^i) H B) = 0, we must have 

'^BiyuVi) gyAy2) fiiu)gu{yi) = crxo-xo- - a.e. , 
which imphes from the hypothesis (fT4l) on g 

li?(l/i,y2) = o- X o- X cr - a.e. . 
However this imphes /i = a— a.e., which is a contradiction. □ 

Proposition 6.3. There is a unique q.s.d. associated with the exponential 
decay rate 6 = X — b. 

Proof. Let u and u' be two different q.s.d. with the exponential decay rate 
9. We can write the Lebesgue decomposition 

^' = /^ + e 

with / a nonnegative measurable function and ^ a singular measure with 
respect to z/. Assume ^ 7^ 0. Applying i?J we get 

„' = fu + ^ = Rl{fu) + Rl{^). 

If / is bounded, since u is a. q.s.d., we have R\{fi') ^ v. In the general 
case, the same result holds by approximating / by an increasing sequence of 
nonnegative functions. Therefore, we must have 

Integrating the function (pi as before, we conclude that R\{fv) = fu, and 
therefore -RJ(0 = C 

Then, we have two q.s.d. u and ^ with exponential decay rate 6 = X — b, 
which are mutually singular. We claim that this is excluded by Lemma [6. 2[ 
Indeed let i? be a measurable subset such that C,{B) = u^B'^) = 0. Then 
i^{B n ^(11)) = iy{A(^ii)) > 0. Since « jj, we get fi{B) > ij,{B fl ^(n)) > 0. 
From Proposition 16.21 we deduce ^{B n v4(ii)) > which is a contradiction. 
Namely ^ = and we conclude that z/' = fu. Let us now show that / = 1, 
which will yield v' = v and so will conclude the uniqueness result. 

Recall the following notation on the ordered lattice of measures ^f{A~^): 
\v\ = f "•" + (— f)^ with f = max(f , 0). Since the linear operator i?J is positive 
it holds < -Rl|f I, that is for all positive and measurable functions (p, 
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it holds / ipd\Rl{v)\ < J (pdRl\v\. Moreover, when there exists a couple 
of sets Ai, A2 such that RI{v){Ai) > > RI{v){A2) and Rl\v\{Ai) > 0, 
i?J 1^1(742) > 0, this inequality becomes strict and we put |-Rj(f)| < -Rj|f|. 
This means that 

V^>0, /^dii;M< 00 where j ^ d\R\(v)\ < j ^ dR'M . 



Assume that u{f 7^ 1) > 0, which implies 7^ 1) > 0. Thus, the sets 
Ai = {f < 1} and = {/ > 1} fulfill the requirements for the signed 
measure v = v — v' . Then, by using that v and v' are R\ invariant, we find 

\v - u'\ = |4M - I < 4(1^ - ^'D ■ 

Since Lfi is Rt invariant and J Lpid\u — u'\ < 00, we find 

y y9i c/[z/ — z/'[ < J dRl{\i^ — = J d\v — v'\ , 

which is a contradiction since v^i > 1. The result is shown. □ 

The proof of Theorem 11.31 is now complete, it follows from Propositions 16.11 
and l6.3[ 
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